Curie temperature of Kondo lattice films with finite itinerant charge carrier density 
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We present a model study of ferromagnetic films consisting of free Bloch electrons coupled to local- 
ized moments (Kondo lattice films). By mapping the local interaction onto an effective Heisenberg 
Hamiltonian we obtain temperature and carrier density dependent exchange integrals mediating the 
interaction between local moments via the conduction electrons. The non-perturbative approach 
recovers analytically the weak-coupling RKKY interaction and yields convincing numerical results 
in the strong coupling (double exchange) regime. The Curie temperature is calculated for various 
coupling strengths, band fillings, and numbers of layers. The results are compared with total en- 
ergy calculations. We discuss the influence of charge transfer between film layers and of anisotropy 
on the Curie temperature. The model we investigate is considered relevant for the understanding 
of the basic magnetic properties of manganites, diluted magnetic semiconductors, and rare earth 
substances as, e.g., Gadolinium. 
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I. INTRODUCTION 

The Kondo lattice model (KLM) is one of the proto- 
type models in solid state physics whenever the coupling 
of itinerant charge carrier spins with immobile localized 
magnetic moments has to be considered. There are vari- 
ous quite different classes of materials for which the KLM 
is used to describe the electronic and magnetic properties 
at least in principle. 

In its simplest single-band version the KLM Hamilto- 
nian reads 



(1) 



(t) 



The first term describes non-interacting electrons (c\ 
annihilates (creates) an electron of spin a =],[ at lattice 
site i) with hopping integrals Ty. er^ is the electron spin 
and Si the localized spin operator, and both are exchange 
coupled by J. 

Depending on the sign of J one is confronted with quite 
different physics. Numerous studies have been carried to 
understand the physics of heavy fermion systems. Here 
the coupling between electrons and localized spins favors 
an antiparallel configuration (J < 0) leading to a com- 
petition between RKKY interaction and spin screening 
effects. 

In this work we want to focus on a favored parallel spin 
alignment (J > 0, ferromagnetic Kondo lattice model, 
FKLM). The FKLM was introduced to model the mag- 
netism in some manganese compounds ("manganites"). 
In the framework of a simple two-site model this kind of 
interaction was found to result from the so-called double 
exchange lending its name to the lattice version in the 
strong coupling regime ~£ 

In the manganites the 5 Mn d-shells are split by the 
crystal field into three degenerate t2g-orbitals forming lo- 
calized spins of S = I which interact via Hund's rule 
with itinerant electrons stemming from the remaining 
two degenerate e g -orbitals. Due to the Jahn- Teller effect 
and the strong Coulomb interaction among the conduc- 
tion electrons the manganites exhibit a complex phase 



diagram, the most prominent phenomenon of which is 
colossal magnetoresistance.^ Although for a more quan- 
titative description of manganite materials additional in- 
gredients like orbital physics, electron-polaron interac- 
tion, and a strong Hubbard-like Coulomb interaction are 
needed, numerous research work on the subject has been 
carried out using the simple one-orbital model (fTj). 3 

More recently attention was directed to the diluted 
magnetic semiconductors (DMS) that are thought to 
have the potential for promising technical applications 
in microelectronics i 5 -* 6 - Substituting transition metal im- 
purities like Mn in a semiconducting host (often III/V 
like GaAs) localized magnetic moments are introduced 
with the consequence of ferromagnetic ordering of the 
randomly distributed maganese cation spins that inter- 
act via their (ferromagnetic, intermediate) coupling to 
valence and impurity band holes. 

Another important field of application of the FKLM is 
the theoretical description of rare earth substances like 
Eu and Gd and their compounds. Here it initially served 
to explain the famous red-shift of the absorption edge 
of the optical 4/ — 5d transition in the ferromagnetic 
semiconductor EuO. 7 Combined LDA/many-body calcu- 
lations based on the FKLM yielded realistic values of the 
Curie temperature for Gd^ and predicted surface states in 
EuO-films 9 . Doping Gd into EuO makes for an attrac- 
tive means to tune the charge carrier density and thus 
the carrier-induced coupling among the magnetic ionsJ^ 
Having the same benefit, Gd-doped GaN was reported to 
yield high Tc ferromagnetism above 300 K^i 

In all of the above-mentioned substances, especially 
regarding possible technical applications in the future, 
it is important to understand and make use of the spe- 
cial physical properties due to reduced dimensionality. 
Several recent works focussed on the role of magnetic 
anisotropy effects in thin films of manganites 12 or of DMS 
films like (Ga,Mn)AsH. Gd films are investigated due to 
their potential employment in spintronics as well as moti- 
vated by the "evergreen" issue of a possibly enhanced sur- 
face Curie temperature . 14 ! 15 ' 16 It is furthermore a well- 
known fact that the Curie temperature strongly depends 
on the film thickness. Due to the broad relevance of the 
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FKLM in current research on magnetism and electronic 
correlations and due to the importance of film structures 
it is worthwhile investigating the dependence of the fer- 
romagnetic transition temperature on the model param- 
eters. 



From the theoretical point of view our treatment of 
FKLM films that is divided up into two parts. First we 
solve the problem for the electronic self-energy within an 
equation of motion approach that is based on a decou- 
pling scheme fulfilling non-trivial limiting cases. In the 
second step we map the interaction onto an effective lo- 
calized spin Hamiltonian by integrating out the electronic 
degrees of freedom. This results in effective exchange in- 
tegrals which incorporate the interaction of the localized 
moments via the conduction electrons beyond conven- 
tional RKKY theory. They depend on temperature and 
carrier density. In sum we end up with a self-consistent 
theory which does not require the assumption of classical 
core spins and which is employable at all temperatures. 



Results on the magnetization and critical temperature 
of ferromagnetic KLM model films with finite band oc- 
cupation have already been reported in Ref. fl7l . How- 
ever, in the present work we make use of a more sophis- 
ticated approach for the fermionic subsystem and carry 
out a wider and considerably more detailed investiga- 
tion of the Curie temperature. A similar system is also 
treated in Ref. [H based on a method guided by the 
same spirit. However, the theory of the electronic sub- 
system there is on a mean-field level only and just nearest 
neighbor Heisenberg interaction is considered, in spite of 
the underlying RKKY-like exchange being long ranged. 
What's more the effective exchange integrals derived in 
Ref. [H in terms of the susceptibility of free electrons are 
proportional to J 2 and thus fail to give the appropriate 
strong-coupling behavior where the Curie temperature is 
independent of J. As far as the numerical evaluation is 
concerned the thickness dependence of the Curie temper- 
ature is analysed. 



This paper is organized as follows. In the next section 
we briefly summarize our method to solve the KLM for 
the electronic self-energy. Then we derive an effective 
Heisenberg Hamiltonian which results from integrating 
out the fermionic degrees of freedom. In the numeri- 
cal results section we discuss the temperature dependent 
quasiparticle excitation spectrum and charge transfer be- 
fore we present an extensive analysis of the Curie temper- 
ature and its dependence on the various model paramters 
like, e.g., the magnetic anisotropy strength. 



II. ELECTRONIC SELF-ENERGY 

In second-quantized notation and introducing Greek 
indices to number the layers the Hamiltonian (TT]) becomes 

H = H a + Hj (2) 

= ^i T ?/ ~ » S if)4 a <T C 30<r ( 3 ) 

ija 

—J \ ZaSftyUiotcr + S ia c\ a(T Ci a -a^j • (4) 

S£~ = Sf a ± iSf a , | (|) = + (-), z T = 1, z x = -1, and 
^ia ' Z are ^ ne cartesian components of the localized spin 
operator Si a . rii a(T is the particle density operator. 

In the following we apply an equation of motion (EOM) 
approach based on a decoupling scheme already applied 
to semiconducting Kondo lattice films (EuO 9 and EuS 19 ) 
and to bulk Gd ? in combined LDA and many-body 
calculations (moment conserving decoupling approxima- 
tion, MCDA). For details of the method we refer the 
reader to corresponding model studiesj 20 i 21 ' 22 After writ- 
ing down the EOM of the one-particle Green function 
G ifa( E ) = (( c ia<y;c f jfjcr )) and the generated higher Green 
functions a decoupling is performed that is guided by 
some non-trivial cases, as, e.g., the ferromagnetic semi- 
conductor at T=0. In the k-space defined within the film 
plane and using matrix notation one can write: 

G ka (E) = [EI - e(k) - VuAET 1 ■ (5) 

e(k) is the hopping matrix. A local approximation of the 
self-energy £°^(i£) — > SijS a fjY^ acr (E) is performed. This 
corresponds to neglecting magnon energies in the three- 
dimensional case 2 ^ and is justified by the fact that these 
energies are usually orders of magnitude smaller than the 
exchange coupling J and the electron bandwidth W. 

From (O one immediately obtains the one-particle lo- 
cal density of states (LDOS): 

Paa (E) = - \lmG™ (E + z0+ - fi) . (6) 

It should be mentioned that in general the MCDA ap- 
proach does not capture Kondo scaling and especially 
violates the Luttinger theorem lm'E er (E = cf) —> 0. 2 - 
Although this might well be due to our approximative 
scheme questions if and in which parameter regimes the 
ferromagnetic Kondo lattice model is a Fermi liquid still 
remain open. Nevertheless we would consider it worth- 
while to investigate the low temperature Kondo physics 
within the framework of our theory. In the present work, 
however, we rather want to focus on the magnetic tran- 
sition temperature. 

The self-energy matrix X CT (i?) depends on various ex- 
pectation values of pure fermionic, mixed fcrmionic-spin, 
and pure localized spin character: 

K = F((n aa ), (5- CT cLc Q - CT ), (S z a n aa ), (7) 

(S z a ), ((S z a ) 2 ), <(^) 3 }, (s^s-)) . 
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The first two types can be calculated self-consistently 
within the MCDA using the corresponding Green func- 
tions and the spectral theorem 2 ^ The localized spin 
correlation functions like, to lowest order, the layer- 
dependent magnetization (S 1 ^) need further considera- 
tion. In order to evaluate these we map the interaction 
term (|4]) onto an effective coupling between the localized 
spins, which is carried out in the next section. 



conventional RKKY interaction is reproduced. In order 
to include spin scattering effectively we dress the free 
propagator in (fT2| and replace it by the full propagator 
Gk+qa/3er(-E)- Within this modified RKKY approxima- 
tion (MRKKY) the exchange integrals in the effective 
Hamiltonian 



H 



ija/3 



(13) 



III. EFFECTIVE LOCALIZED SPIN 
HAMILTONIAN 

It is well known that using perturbation theory (j4]) 
leads to the so-called RKKY interaction^^ This second 
order interaction between localized spins is long ranged 
and yields an oscillatory behavior of the exchange cou- 
pling which is mediated by uncorrelated electrons. It 
appears highly desirable to have such a mapping with- 
out being restricted to the weak coupling regime. This 
requires finding an effective Hamiltonian 



H 



eff _ 



(Hj) 
J 

~2N 



e" iqR? (S, Q • a) 



aa' \ L -k+qao- C kao- 



(8) 

\(c) 



are given by 



+ 00 



k<J -oo 

(14) 

where / is the Fermi function and the sum is over the 
first Brillouin zone. 

Within the Tyablikov approximatio n 27 ! 28 the imagi- 
nary part of the transversal layer-diagonal spin Green 
function (spectral density) can be written as 

- "Mfe k)a » = 2(S' a ) J2 €^ (E ~ £r(k)I15) 



where N is the number of lattice sites in one layer and 
the superscript c formally indicates that the averaging is 
performed in the fermionic subspace only. The method 
we present in this section has been applied before to bulk 
Kondo lattice systems^^ We generalize it into a matrix 
formulation which allows us to consider film structures, 
but is also relevant whenever sublattice decompositions 
are necessary, for instance when dealing with antiferro- 
magnetic configurations. 

In order to obtain the expectation value in ((8|) we in- 
troduce the modified Green function 



(E) = ((c kQ0 



-k+q/3 CT 



(c) 



(9) 



The EOM for G reads 

; {E8 ia - G& + ^(E) = <5 q0 <W'<W (10) 
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J 

'2N 



J- e *( k -PW ( SiQ • a) a , a „ G^ W (E) . (11) 

ipa" 

In principle, this equation can be solved iteratively. Ap- 
plying the spectral theorem then yields the expectation 
value in ([8]). This correlation function is an operator 
since we are working in the fermionic subspace, resulting 
in an effective interaction among the localized spins 
only. 

We have to find a manageable approximation for G on 
the right hand side of (fT0| . It can easily be shown that 
in the non-interacting limit, 

GpM+qapiE) — > ^tr" crSpM+qG^l^^E) , (12) 



where the spectral weights rj^ 1 and energy poles -E 7 (k) 
have to be evaluated numerically. For a monolayer we 
can drop all Greek indices, ryk = 1, and we get the well- 
known spin wave energies 



£(k) = 2(^)(J(0)-J(k)) 



(16) 



Applying the Callen metho d 27 i 29 to superlattices the 
layer-dependent magnetization, 

/qz , _ (i + y a + s)^ a s+1 + (s - ya )(i + ^ a ) 25+1 n 71 

[ a> (1 + v*) 2S+1 < S+1 1 ' 

and other higher order spin correlation functions can be 
obtained using the Bosc-likc distribution function 



1 



k.7 



i 



(18) 



([5)) . (fT4|) . and (fT7|) represent a self-consistent system of 
equations that can be solved for the one-particle Green 
function matrix Gfc a (E) and the magnetization (S^). 

Before proceeding to the numerical evaluation we have 
to reconsider our effective Hamiltonian (fT3"|) . It is known 
that for low-dimensional systems anisotropies can be- 
come very important and even a necessary condition for 
magnetic ordering at finite temperature i 30 ! 31 We there- 
fore include a single-ion anisotropy term 



(19) 
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The physical background of this magnetic anisotropy is 
spin-orbit coupling, which usually is some order of mag- 
nitudes smaller than the exchange coupling between lo- 
calized spins. A positive K% favors an out-of-plane easy 
axis, i.e. perpendicular to the film plane. 

We treat this term in the Anderson-Callcn 
approximation^ and decouple the higher Green 
function generated by (|19|) in the following manner: 



[Sf,H A ]_;ST)) 



+ & facial Sjp)) 



$ Q = 2(5*) ( 1 - — (5(5 + 1) - <(5*) 2 }) 



(20) 



(21) 



We can drop the site index i due to lateral translational 
invariance. The single-ion anisotropy acts as an effective 
field <& Q coupled to 5„ . Note that this decoupling is valid 
only if the magnetization is parallel to the z-direction. 
In order to ensure this a rotation of the coordinate sys- 
tem might be required before the decoupling, for instance 
when an arbitrary oriented magnetic field is present^ 

For bulk or a monolayer (a = 1) it is possible to derive 
a simple formula for the Curie temperature Tq. Expand- 
ing (fl7|) and (fT8)) in the vicinity of Tc one easily finds 




0.5 

E-\L [eV] 

FIG. 1: One-particle density of states of a 5-layer film at dif- 
ferent temperatures. Parameters: n — 0.2, J — 0.2 eV, S = 
7/2. Critical temperature Tc = 174 K. Charge neutral calcu- 
lation. 



A. Electronic properties 



1. One-particle excitation spectrum 



S(S + 1)N 



E 



tf 27 + 2(J(0)-J(q)) 



(22) 



T=T C 



where 



lim 



$(T) 2(25-1) 



T^Tc (5 2 > 



35 



(23) 



is a constant and depends on the specific decoupling of 
the anisotropy Green function (f2"D")) j^ 

For a > 1 the spin wave energies are no longer known 
in analytical form as in (|16[) . One can then, however, use 
(|17p to determine the critical temperature. 



IV. NUMERICAL RESULTS 

All numerical results have been obtained for simple 
cubic (sc) (100) films. The nearest-neighbor hopping 
integral t in the tight binding approximation was cho- 
sen according to a bulk bandwidth W — 1 eV, i.e. 
t = —0.083 eV. The value of the exchange coupling is 
assumed to be homogeneous and isotropic for all sites 
and layers of the films, as already implied by l[4"]). Fur- 
thermore for the anisotropy constant we set K% = K2 
for all layers. The magnitude of the localized spin is 
5 = 7/2. To keep notation simple we write ni a>a ) for 
the respective expectation values (nr aj(7 \}. n denotes the 
average electron density. 



We start with the discussion of the local one-particle 
density of states. Fig. Q] shows the LDOS at different 
temperatures between T = and Tq. The Curie tem- 
perature has been self-consistently computed. 

The general picture for a layered system is the same 
as for the bulk Kondo lattice models At low tempera- 
ture the spin-up spectrum basically consists of the rigidly 
shifted non-interacting electron energy band. At T = 
and in the insulating limit (band occupation n = 0) this 
is a rigorous result^ For the spin-down LDOS the sit- 
uation is different. Besides a band at higher energies 
corresponding to an antiparallcl coupling of conduction 
electrons and localized spins there is a scattering part in 
the energy range of the f-band reflecting the fact that a 
^-electron can flip its spin by creating a magnon in the 
localized spin subsystem. With increasing temperature 
spin symmetry in the spectrum is gradually established. 
This is accompanied by a reduction of the bandwidth due 
to reduced effective hopping by spin scattering. For in- 
termediate coupling strengths J ss W as in Fig. [1] this 
leads to a temperature-induced opening of a gap in the 
excitation spectrum. From the zero bandwidth limit one 
can learn that the distance between the two correlated 
bands roughly scales as ~ J(5 + \). 

Apart from these facts already known from bulk re- 
sults one also observes typical features of reduced dimen- 
sionality. The number of van-Hove singularities which 
are most pronounced at low temperatures is indicative of 
the number of layers and therefore for the finiteness of 
the film system. A second important observation is the 
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FIG. 2: Charge transfer in a 5-layer film at T = 0, (S z ) = S 
(ferromagnetic saturation) . Whereas the carrier density in the 
surface layer is always reduced in the non-interacting system 
there is an occupation above average for large enough n in 
the correlated system. 
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FIG. 3: Charge transfer in a 5-layer film at T = 0, (S z ) = 0. 
The horizontal line at n = 0.75 indicates the total average 
occupation number. The choice of the energy scale refers to 
the free simple cubic bulk bandwidth W = 12t. 



reduced effective bandwidth of the surface compared to 
the center layer LDOS. This can be directly traced back 
to the lower variance of the surface LDOS: in the non- 
interacting limit A 2 p = qt and q sur face/qbuik = 5/6 for 
a sc(100) geometry. The effect is thus not correlation in- 
duced and more pronounced for "open" film geometries. 
The band edges of the surface and center LDOS, however, 
are the same. 



2. Charge transfer 

The broken translational symmetry is furthermore re- 
sponsible for the occurence of charge transfer, i.e. the 
deviation 5 a n of the layer dependent band filling from 
the average occupation number: 

S a n = n a - n = ^(n Qcr - — ^n 7cr ) . (24) 

cr 7 

Nl is the number of layers. Charge transfer is already 
present in the free electron film system and is due to 
the smaller effective bandwidth of the surface layer (oc- 
cupation n s ). From this follows directly that in order to 
ensure thermodynamic equilibrium, i.e. a common chem- 
ical potential /i, one has n s < n below and n s > n above 
half-filling (n = 1). 

Switching on the coupling of the electron spins to 
the localized spins the question arises how this behav- 
ior changes. Fig. [2] shows the difference of the surface 
and the average occupation number as a function of n at 
T = and ferromagnetic saturation. Whereas the elec- 
tron density in the surface layer n a —i is always below the 
average density n in the non-interacting limit (J = 0) one 
also finds 8 a n > at n < 1 for J > 0. The explanation 
follows the same reasoning as above for the free system: 



a reduced or enhanced occupation in the surface layer is 
linked to the chemical potential being in the lower or up- 
per half of the LDOS, respectively. As there is a (mainly 
J.-) transfer of spectral weight to higher energies for fi- 
nite J the chemical potential crosses the maximum of 
the lower band at values n < 1. One finds the quali- 
tatively same and quantitatively similar behavior in the 
spin disordered phase where (S°) = 0. 

A second point is the dependence of the charge trans- 
fer on the coupling strength J. The interaction term Q 
alone energetically favors single occupation of a lattice 
site for n < 1. This can be easily seen in the zero band- 
width limit. 25 ' 35 It can thus be expected that a finite J 
suppresses charge transfer. This is indeed the case as Fig. 
[3] shows where the electron density for the different layers 
of a 5-layer film is plotted as a function of the exchange 
coupling. The change of sign of 6 a n is again caused by 
the same mechanism as pointed out above. Given a band 
occupation of n < 1 the chemical potential moves from 
the lower half of the Bloch bands (J = 0) to the upper 
half of the lower quasiparticle subband (J > 0). 

We observe that even J — ► oo does not lead to com- 
plete charge neutrality and that at large J the higher- 
occupation of the surface layer goes with an almost ho- 
mogeneous charge distribution in the inner layers of the 
film. It is also interesting to note that at J w 0.1 there is 
a quasi-homogenous charge distribution which coincides 
with the change of sign of 5 a n. 

One can summarize that the effect of charge transfer 
is highest at low couplings J and low average band oc- 
cupations n, where it can amount to almost 30% (Fig. 
[5] at n = 0.1). It appears to be of minor importance in 
the strong coupling regime, especially at higher charge 
carrier densities. 

All the findings we have discussed in this section are 
qualitatively the same for other film thicknesses. As al- 
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FIG. 4: Curie temperature for various numbers of layers. Pa- 
rameters: n = 0.2, K2 = 10 -6 eV. There is no significant in- 
fluence of charge transfer on the critical temperature. Inset: 
Curie temperature of a monolayer for a low band occupation. 



FIG. 6: Total, kinetic, and potential energy differences AU = 
Upm(T=0) — Ufm(T=0) between the paramagnetic and the 
ferromagnetic state as a function of the exchange coupling 
J for a monolayer. The negative values of A U at small J for 
n = 0.6 indicate that ferromagnetism is unstable. 



ready mentioned charge transfer effects are more pro- 
nounced in "open" surface geometries with a stronger 
relative reduction in the number of nearest neighbors of 
a surface atom. However, the order of magnitude of the 
effect is the same. 

Of course in real systems, Coulomb interaction pre- 
vents charge separation on a macroscopic scale. For 
the present purpose it is sufficient to discuss how charge 
transfer affects the ferromagnetic transition temperature 
of Kondo lattice films. To this end we will compare 
results based on calculations with and without charge 
transfer. Enforcing charge neutrality means we evalu- 
ated the layer-dependent centers of gravity of the Bloch 
bands self-consistently such that n a = n for all a. Where 
not otherwise stated the following results were obtained 
for charge neutral KLM films. 
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J/W 

FIG. 5: Curie temperature for various numbers of layers. Pa- 
rameters: n = 0.6, K2 = 10 -6 eV. 



B. Curie temperature 

1. Whole J -range 

Our non-perturbative theory allows for an evaluation 
of KLM films at all coupling strengths J. In what follows 
the anisotropy constant is much smaller than the other 
energy scales in our model, K% << W, J. The influence of 
the anisotropy on the Curie temperature will be discussed 
in the last subsection in more detail. 

Figs. [4] and [5] display the ferromagnetic critical tem- 
perature as a function of the exchange coupling J for 
various film thicknesses. The inset in Fig. 2] shows Tq 
of a monolayer for a small band occupation n = 0.01. 
Here one recognizes the typical Tc ~ J 2 -behavior of the 
perturbational RKKY interaction for small J. Increas- 
ing the electron density induces a critical interaction J c 
below which ferromagnetism is not stable. For a given 
number of layers J c increases with n. It has been shown 
before analytically that there is a minimum critical in- 
teraction strength for stable FM in the KLM in infinite 
dimensions^ 

To discuss this point further we consider the internal 
energy per site U: 



U=(H) = (H Q ) + (H. ; ) = Uun + U pot 
It can easily be calculated with the LDOS ([6]): 



U = — 



Nr ^ 



+ 00 



dEf(E)Ep^(E) 



(25) 



(26) 



We start from F = U — TS where F is the free en- 
ergy and S is the entropy. By comparing the differ- 
ence of the internal energy at T=0 between the ferro- 
magnetic (FM) and paramagnetic (PM) state, AU = 
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FIG. 7: Total, kinetic, and potential energy differences be- 
tween the paramagnetic and the ferromagnetic state AU = 
Upm(T—Tc) — Ufm(T—0) as a function of the band occu- 
pation for a monolayer in the strong coupling regime. The 
Curie temperature obtained with the MRKKY theory is also 
included for comparison. Inset: MRKKY Curie temperature 
for various numbers of layers. Parameters: J = 1 eV, K2 = 
1(T 6 eV. 
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FIG. 8: Total, kinetic, and potential energy differences be- 
tween the paramagnetic and the ferromagnetic state AU = 
Upm(T—Tc) — I7fm(T=0) as a function of the Curie tem- 
perature in the strong coupling regime. The film thickness 
is an implicit parameter and increases from left (monolayer) 
to right (Nl = 15). Parameters: n — 0.5, J = 1 eV, K2 = 
1(T 6 eV. 



Upm(T=0) — Ufm(T=0), we can evaluate the stability of 
the ferromagnetic state against paramagnetism. In Fig. 
O the kinetic, potential, and total energy differences are 
shown for a monolayer at T=0. The magnetic stability 
for low to moderate J is governed by a complex interplay 
between [/kin and U po t- Whereas at n = 0.2 it is the po- 
tential energy which favors the (ferro)magnetic state it is 
both energies which make the ferromagnetic state unsta- 
ble for higher n, leading to a critical interaction strength 
as in the self-consistent calculation based on the effective 
Heisenberg model before. The fact that in terms of the 
ground state energy there is no J c at n = 0.2 contrary to 
the results in Fig. 2] suggests that other magnetic corre- 
lations such as non-commensurate or antiferromagnetic 
become important, reducing further the parameter range 
of FM stability. 

Furthermore it is remarkable that the critical interac- 
tion strength is smallest for a monolayer whereas films 
with a > 1 behave rather bulk-like as far as J c is con- 
cerned. Taking into consideration the energy scales of 
the KLM one would expect this tendency. Quite gener- 
ally the physics should be governed by the ratio J/W. 
Hence a smaller bandwidth due to a reduced number of 
layers implies a smaller J around which the intermediate 
coupling regime is located. In addition the reduced bulk 
FM region might indicate favored antiferromagnetic con- 
figurations which are not existent in the 2D case for mere 
geometrical reasons. 

For comparison we performed calculations where 
charge transfer was permitted. Allowing charge transfer 
means that we used equal center of gravities of the non- 
interacting local density of states Tff = To for all layers. 
We did not find any noteworthy influence of charge trans- 



fer on Tq compared to the charge neutral case as Fig. 0] 
demonstrates for a 5-layer film.— 

2. Double exchange regime 

In the strong coupling (double exchange) regime 
(JS » W) one observes that Tc runs into saturation 
for large enough J at all electron densities and numbers 
of layers. Once more we will contrast the results ob- 
tained by the MRRKY method with energy considera- 
tions. Whereas before we were interested in the stability 
of ferromagnetism at T=0 we want to be more quanti- 
tative in the following comparison of the Curie tempera- 
tures. Hence we evaluate the internal energy of the para- 
magnetic state at Tc as the energy scale for the critical 
temperature is set by the change in total energy when 
going from the ferromagnetic to the paramagnetic state 
T c ~ AU = AU(T C ) = Upm(T=T c ) - U FM {T=0). In 
the double exchange regime Tc is essentially determined 
by the kinetic energy: in the local frame the electron spin 
is oriented parallel to the localized spin and the poten- 
tial energy does not change much between the ferromag- 
netic and the paramagnetic state; for J — > 00 it stays the 
same. On the other hand hopping of electrons strongly 
depends on the magnetic configuration and is somehow 
blocked when the localized spins are disordered. For the 
bulk KLM the relationship between Tc and A [/kin has 
already been analysed ] 37 ' 39 

The inset of Fig. [7] shows the Curie temperature as 
a function of the band occupation n for strong cou- 
pling. We obtain ferromagnetism for a wide range of 
n up close to half-filling and a symmetry approximately 
around quarter- filling for all films. Tc increases with the 
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FIG. 9: Curie temperature as a function of the film thickness. 
Lines are guides to the eye. Parameters: n = 0.2, Ki = 
I(T 6 eV. 



number of layers. The maximum band occupation for 
which ferromagnetism exists does not depend on the film 
thickness, indicating that the FM phase boundary in the 
strong coupling regime is the same in 2D and 3D. Also in 
Fig. [7] a comparison between AC/ and Tc is shown. The 
variation of the critical temperature with n is paralleled 
by a corresponding variation of At/kin- At/ po t varies only 
slowly with n. We find that, apart from a small differ- 
ence in the maximum n yielding ferromagnetism, there is 
even quantitatively a good agreement between the values 
of T c and AC/. 

In Fig. [5] we show AC/ki n and AC/ po t as functions of Tc 
with the number of layers as implicit parameter. One 
recognizes a linear dependence between the ferromag- 
netic transition temperature and the energy change. Tc 
is again dominantly determined by the difference in ki- 
netic energy of the FM and of the PM phase. We point 
out that the almost quantitative agreement of yMRKKY 
and AC/ for the monolayer is due to the particular choice 
of the anisotropy parameter and thus rather coincidental. 
With an increasing number of layers deviations between 
the two quantities appear. However the magnitudes re- 
main similar. 



3. Tc as a function of film, thickness 

We conclude this section by presenting the critical tem- 
perature as a function of the film thickness. As a natural 
consequence of the results and the discussion just pre- 
sented this dependence varies drastically in the different 
coupling regimes as demonstrated in Fig. [5] Whereas for 
large J the Curie temperature decreases when reducing 
the number of layers, reflecting the findings just discussed 
before in Fig. [5J we get a completely different picture for 
weak to intermediate strength of J. At J = 0.16 the fer- 
romagnetic transition temperature is equal for all films. 




0.4 0.6 

J/W 

FIG. 10: Curie temperature of a monolayer at different 
anisotropy strengths for n = 0.6 (from bottom to top: K2 = 
1, 10, 100 fieV). Varying the anisotropy has a considerable 
effect on Tc in the strong coupling regime but only slightly 
changes J c . Inset: Curie temperature as a function of K2 for 
a monolayer in the saturation region (n = 0.2, J = 0.2 eV, 
see also Fig. [4}. 



Reducing further the exchange coupling J yields a finite 
Tc only for very thin films. 



4- Anisotropy 

In the last subsection we want to briefly discuss the ef- 
fect of a single-ion anisotropy on the ferromagnetic tran- 
sition temperature of Kondo lattice films. There exists 
abundant literature on scaling relations between the crit- 
ical temperature and key film parameters as anisotropy 
strength, exchange coupling, and the film thickness, 
mostly for Ising-like or Heisenberg films. As we are deal- 
ing with an effective localized spin model, derived from 
the KLM, these results should also apply in our case. The 
exchange integrals (fT4|) now are, however, T-dependent. 
In the following investigation of the functional depen- 
dence Tc(J,K%) we restrict ourselves to a monolayer 
because in this case the role of anisotropy is most pro- 
nounced and the extension to multilayers is straightfor- 
ward. 

The effect of the anisotropy strength on the ferromag- 
netic transition temperature over the whole range of the 
intra-atomic exchange coupling J is shown in Fig. 1101 
One can distinguish two regimes: As can be seen from 
the inset (representing the saturation regime) the well 
known logarithmic dependence of the Curie temperature 
on the anisotropy strength K2 is also valid for the Kondo 
lattice model. Already a very small anisotropy is suffi- 
cient to obtain a Tc which is of the order of magnitude 
of Curie temperatures realized for much larger, more re- 
alistic values of K%. 

Evaluating the formula for the Curie temperature (|22p 
by solving the corresponding q-integral for the predomi- 
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nantly contributing small wave vectors one gets 

^"( 1+c %r)- (27) 

c is a constant w 1 and J is the T— dependent effective 
coupling of one localized spin to all others, i.e. it de- 
termines the order of magnitude of the magnon energies, 
which are typically much larger than the anisotropy en- 
ergy, J >> K 2- An analogous formula for the Heisenberg 
model with constant exchange coupling is discussed, e.g., 
in Ref. Hcl In our case the T-dependence of the exchange 
integrals IT4|) in the paramagnetic regime is evidently too 
weak to cause significant deviations from the logarithmic 
dependence, not too surprisingly because this tempera- 
ture dependence is a mere Fermi softening effect. 

Whereas in the saturation region the scaling of Tc with 
K 2 is determined by Eq. (f27|) as just discussed, in the 
critical region FM becomes unstable and a higher K 2 only 
slightly reduces J c . There Tc — > and thus the effective 
coupling J — > 0, see Eq. (|2"2")l where the effective coupling 
corresponds to the J-exchange term. The critical region 
has a width of AJ s=s 10 _1 eV within which the effective 
coupling typically changes by AJ w 10~ 3 eV. Thus a 
AK 2 ~ 1CT 4 eV should shift J c by w 1CT 2 eV, which is 
in accordance with the numerical findings in Fig. 1101 

V. CONCLUSIONS 

We have investigated ferromagnetic Kondo lattice films 
with a finite band occupation. Our self-consistent theory 
is made up by an electronic and a magnetic part. The one 
particle Green function is obtained within an equation of 
motion approach. Local moment correlation functions 
are calculated by an effective Heisenberg Hamiltonian 
which results from integrating out the electronic degrees 
of freedom. 

We presented results for the temperature dependent 
one particle exciation spectrum and discussed charge 
transfer between film layers. Charge transfer occurs in 
all coupling regimes and, at strong coupling, leads to a 
pronounced deviation of the surface layer electron den- 
sity from the rather homogeneous band occupation of the 
inner layers. The reduced dimensionality is also responsi- 
ble for a reduced magnetization at the surface. However, 
Tq is the same for all layers. 



Curie temperatures were discussed for a variety of pa- 
rameters like the exchange coupling J, the band occu- 
pation n, and for various number of layers. For large 
enough band occupation we obtain a critical J c below 
which ferromagnetism is not stable. In the strong cou- 
pling regime Tc saturates and increases with the number 
of layers. We have not found any significant influence of 
charge transfer on the Curie temperature of KLM films. 

In addition to the MRKKY self-consistent X^-results 
energy calculations of the paramagnetic and ferromag- 
netic state were used to discuss both the critical inter- 
action strength in the intermediate coupling regime and 
the crucial role the kinetic energy plays in the strong cou- 
pling regime. Depending on the values of the exchange 
coupling and the band occupation the dependence of the 
Curie temperature on the film thickness exhibits quite 
a different behavior. We also discussed the effects of 
a single-ion anisotropy on the ferromagnetic transition 
temperature. Whereas it can significantly elevate Tc in 
the strong coupling regime it hardly affects the critical 
interaction strength J c . 

In this work we took the rather classical spin quantum 
number 5=7/2. This would apply for rare earth systems 
as, e.g., Gd. Systems of current interest as DMS (5=5/2) 
possess a somewhat high 5, and even for the manganites 
with 5=3/2 classical spins are believed to be a reasonable 
approximation at least at low temperatures^ Neverthe- 
less a more profound investigation of the dependence of 
the magnetic properties on the spin quantum number at 
higher temperatures appears to be highly valuable. 

A worthwhile extension of the present work would also 
be to use other lattice structures and to take additional 
magnetic configurations like antiferromagnetism into ac- 
count. A combined bandstructure many-body calcula- 
tion for Gd could contribute to the much debated issues 
of an enhanced surface Curie temperature and electronic 
surface states in this substance. More on the model level, 
making the exchange and anisotropy parameters layer 
dependent or analysing the influence of non-magnetic 
cap layers on the magnetic stability of KLM films would 
be interesting tasks. Furthermore an application to the 
problem of interlayer exchange coupling appears to be 
rewarding because the temperature dependent effective 
exchange integrals give direct access to the coupling be- 
tween the layers in a film system. This investigation is 
planned for future work. 
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